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PACS 72 . 25 . Rb - Spin relaxation and scattering 

PACS 71.70.Ej - Spin-orbit coupling, Zeeman and Stark splitting, Jahn- Teller effect 
PACS 03 . 65 . Vf - Phases: geometric; dynamic or topological 

Abstract. - We study zero temperature spin dynamics of a particle confined to a ring in presence 
of spin orbit coupling and Ohmic electromagnetic fiuctuations. We show that the dynamics of 
the angular position 6{t) are decoupled from the spin dynamics and that the latter is mapped to 
certain correlations of a spinless particle. We find that the spin correlations in the z direction 
(perpendicular to the ring) are finite at long times, i.e. do not dephase. The parallel (in plane) 
components for spin i do not dephase at weak dissipation but they probably decay as a power 
law with time at strong dissipation. 




Introduction. — Due to recent advances in semicon- 



ductor technology, it became possible to isolate and nia- 
\^ nipulate spins of individual electrons [HE]- For efficient 
T— I . spin manipulation, however, slow spin decay is needed. 
' Spin decay in mesoscpopic devices is generated by two ma- 
, jor sources: Hyperfine interaction with nuclear spins 3' is 
' responsible for spin decay in most materials. However, 
: spin-orbit (SO) coupling can also induce spin relaxation, 
^ and under certain conditions, phonon- [1] or electromag- 
^ , netic field-induced SO relaxation J\ can dominate the de- 
^ ' cay [5 1 6) . As shown in Ref. j7| two-photon (or two-phonon) 
. processes lead to geometrical spin relaxation even in the 
absence of external field and, as pointed out recently, this 
mechanism can become even dominant in hole-doped sys- 
tems 015]. 

Here we make an attempt to understand, whether the 
above-mentioned geometrical spin relaxation can survive 
even at T = temperature. Although Ohmic electromag- 
netic fluctuations were found to lead to a vanishing spin 
relaxation rate at T = [7 , the results of Ref. [7] are 
not conclusive, since they allow for non- exponential leloK- 
ation, common in Ohmic systems. To address this issue 
more rigorously, we consider a ring geometry. Studying 
a ring is, however, not of pure theoretical interest; high 
quality semiconductor rings jlOUll] can in fact also be 



used as quantum spin qubits [TT] , and usefulness of these 
devices depends on spin dephasing, a topic under active 
experimental study [^[6ll8l[T2lfT3] . 

There are two types of spin-orbit coupling in two- 
dimensional electron systems: the Rashba interaction in- 
duced, e.g., by an electric field perpendicular to a two- 
dimensional (2D) layer [2] , and the Dresselhaus coupling 
induced by bulk inversion asymmetry |15) . Our aim is to 
study, how these couplings influence spin coherence for 
an electron confined to a ring, in the presence of Ohmic 
fluctuations. We shall flrst derive the appropriate Hamil- 
tonian for a confined electron, and show that the presence 
of the spin does not influence the orbital motion of the 
confined electron, which is governed exclusively by fluctu- 
ations of the external electric field. The dynamics of the 
spin, on the other hand, is determined by the orbital mo- 
tion of the electron, and has a topological character. We 
find that for weak dissipation the spin does not dephase, 
but certain spin components are reduced by fluctuations. 
For strong dissipation, however, we find that certain com- 
ponents of the spin probably relax even at T = tem- 
perature, due to the disordering of the orbital degrees of 
freedom jl6| . The relaxation we find is, however, not ex- 
ponential but of a power law, typical of confined particles 
at temperature T = pJllTS]. 
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Hamiltonian. — Let us start by projecting a 2D spin- 
orbit Hamiltonian on a ring, a procedure which is not en- 
tirely trivial [T9[|20j . In addition to the kinetic terms, the 
2D Hamiltonian consists of a potential Vq (r) that confines 
the particle to a ring of radius R ± SR, with SR <C R. 
We write the total Hamiltonian in polar coordinates as 
Tio + 'H' where 
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I, S are spin operators, mg is the elec- 



H' = 
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tron mass and Px and Py denote the x and y components 
of the momentum. The ao term is the Rashba coupling 
while Po denotes the Dresselhaus coupling. Labeling the 
radial eigenstates of Ho by \n) and their energies by En, 
our aim is to project H' on the subspace, |0), while keeping 
terms up to order 0{SR), a procedure that involves some 
subtleties. First we rewrite H' = Hi+H2 by introducing 
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As noticed by Meijer et al. [19], for any state ip{r) that 
is radially localized near R, one has {ip\2dr + f IV') — 0. 
Therefore, to first order in the SO coupling, H2 does not 
give a contribution to the projected Hamiltonian. Never- 
theless, as previously overlooked, H2 cannot be ignored: 
localization on a scale SR implies dr ^ l/SR and hence 
2nd order perturbations in H2 do give a contribution, 
- Hi/ En = 0(1), since En - 1/iSR)^. The next or- 
der contributions scale as H^/E"^ = 0{SR), and vanish in 
the limit — )■ 0, similar to all higher order terms in the 
perturbation series. 

Perturbation theory to 2nd order yields therefore the 
projected spin and angle dependent effective Hamiltonian 

^ {0\H2\n) {n\H2\0) 



H 



ring 



(0|Hi|0> 



En — Eq 



0{5R) . (3) 



The sum in Eq. ([3]) can be evaluated analytically by mak- 
ing use of a sum rule [2^11^ , and the second term of Eq. ^ 
simply becomes ime(aoSe — /3oS^)^. 

Introducing the vector h{6) via h = [acosO — 
/3sin6',asin0 — /3cos6'), and with the dimensionless 
Rashba and Dresselhaus couplings defined as a = mRao 
and /? = mRPo, we can finally rewrite our effective ring 
Hamiltonian in the 5R limit as 

t2 



2mgi?2 



[pg + h{e) ■ sy 



(4) 



We remark, in particular, that the term ^ a/3 sin 26 in the 
effective Hamiltonian of Ref. [5D] is exactly canceled by 
the 2nd order terms. As a consequence, Eq. (|4|) possesses 
a conserved "momentum", Q = pg + h{9) ■ S. 
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Fig. 1: Evolution of the spin (coherent state) while the electron 
makes a circle, 60 — > 6o + 2tt. The initial state is rotated around 
an axis N(0o) by an angle Tog. 

Spectrum. — The eigenstates and eigenenergies of ^ 
can be analytically computed for /3 — 0, when the system 
is rotationally invariant and therefore Jz = pg + Sz is also 
conserved. The Hamiltonian can then be written as 



Hri 



2m7W 



[Jz ~ n{9) ■ SVTT^ y 



(5) 



with n{0) = (—hx{9),—hyl6), l)/\/l + a unit vector. 
The energy spectrum and the eigenvalues can then easily 
be found by constructing common eigenstates of the two 
commuting operators, Jz andn(0)-S. For 5 = 1/2, n{9)S 
and Jz have eigenvalues n{9) ■ S = ct/2 and Jz = m + 
(t/2, respectively, with ct = ± and m an integer. The 
spectrum is Cm^a = 2ml r'^ ["^ '^(^ ~ ^ym?)] , and 
the eigenstates are of the form Jg S-=\9) (g) | ± n{9)), with 



I ± n) denoting spin coherent states, defined through the 
usual relation, SI • S|J1) = S\Vt) [24j. The wave functions 
can be explicitly expressed as 



V'm,+ (6') 



( COS f , 

( 



sm 



f ) 



sm ^ , cos ^ 



(6) 



with a defined as a = arctan(a). The states V'im.i are 
related by time reversal, and their energies equal, Em,+ = 
E^rn.-- For a < \/3 the ground state has m — 00 

Dissipation. — Having understood the properties of 
an isolated ring, we now couple the motion of the particle 
to the coordinate ^ of a dissipative environment, i.e., we 
consider the total Hamiltonian as H = Hring + 
Throughout most of this paper we shall assume that 
V{9, ^) describes the coupling to a Caldeira-Leggett (CL) 
environment, appropriate for small rings in an Ohmic 
(metallic) environment!! Then ^ represents the random 
force generated by the environment, V = ^_e*^ -|- ^+e~*^, 
and the T = Fourier transform of the environment cor- 
relations is {£,-£,+ + 5+^_)tj = h"^ 'rj\Lo\, with 77 the dimen- 
sionless friction coefficient. 

The corresponding equations of motion for 9{t) are 



hS 



1 



dgV{9,i) 



^ Here we used the phase convention of Ref. | 24| . This construc- 
tion can be generalized for larger spins. 

^In a dirty metal environment, e.g., one needs the ring's radius 
to be smaller than the mean free path. 
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Hence, as a consequence of the simple form of Hring, 
Eq. (HI, the dynamics of 9 in the dissipative environment 
are not affected by the spin-orbit coupHngs. This decou- 
phng aUows us to describe the 9(t) evolution by a path 
integral, where for each trajectory the spin dynamics fol- 
low from Eq. (g]) 

'^^eh{0)^^ => '^ = hie)xs. (8) 

Viewing as a "time" variable, these dynamics correspond 
to a spin precession around a "time" dependent magnetic 
field h{9). Note that switching to " Schrodinger" picture, 
the spin coherent states have a simple 9 evolution, too. 
Apart from a phase, they evolve as \fl{9)), where n{9) 
is the vector solution of (|8]), i.e. ^ = h{9) x Jl. In 
particular, the vector n{9) can also be shown to satisfy 
this equation. 

In terms of the spin operators, Eq. (jS]) is solved as a sim- 
ple linear mapping Si{9) = Rij{9,9Q)Sj, with Rij{9,9o) a 
rotation matrix. The rotation matrix -R2Tr(^o) = -R(^o + 
2n,9o), corresponding to the particle going once around 
the ring by 2tt, is of special interest. We denote by the 
unit vector N(^o) its axis of rotation and by Tg^ the cor- 
responding rotation angle (see Fig. [T|). In particular, for 
/3 = we find that the angle Fg,, is independent of the 
initial value, F = 2tt(1 — v'T+"c?), and is typically incom- 
mensurate with 27r. 

Mapping to a spinless system. — For a given evo- 
lution, 9q 9, we can obtain the evolution of the spin 
part of the wave function from Eq. ([S]), which is described 
by a unitary operator, [/spin(^, ^o)- Here the Hamilto- 
nian to describe the 9 ("time") evolution of the spin is 
Hg = h{9) ■ S. We proceed to study the case /3 = 0. Then, 
as in the standard NMR rotating field problem, the spinor 
transformation -0' = e*'-^"^"-''^^ V' to the "rotating frame" 
cancels the "time" (9) dependence, and amounts in re- 
placing Hs ->■ h{9o) • S - 5^ = -Vl + a2 n{9o) ■ S. For 
S = 1/2 this leads to the evolution operator, ^7spin(^', 9o) = 
e~*~2~'^'e*^ — 2— nifoj o-_ Using now the expression of 
Q we find that the 9 evolution of the spin states has a 
particularly simple form 

t/spi„(^,0o)0m,±(^o) = e^^-if^"-^) V'™,±(e) , (9) 

where qm ± = 'TT'i|T| \/T+"c? denote the eigenvalues 
of the momentum Q. After a 27r rotation the state ^m,± 
picks up an incommensurate phase, 2TTqm.±- Note that the 
semiclassical evolution involves a similar incommensurate 
angle, F, as discussed below Eq. ^ (see also Fig. 1). 

Making use of the decoupling of orbital and spin de- 
grees of freedom, we can construct a mixed path integral 
formalism (to be detailed in Ref. [23]), where the spin is 
treated in an evolution operator formalism, while the or- 
bital motion of the particle is developed in a path integral 
formalism. The full evolution for a given environment his- 



tory is then obtained as; 

i'm,±{St,t) = y] / d9o / m 

„ "'0 JOo 

e'Spie,0 UspUet + 2^n,9o)^^^±{9o) . (10) 

Note that 9 in this equation is a non-compact variable, 
and an additional integration over the environment con- 
figurations has to be carried out in the end. Importantly, 
the action Sp{9,^) = j^[\mer'^9^ — V{9,C)\ describes a 
particle on the ring in the presence of dissipation for a 
given environment history, and is independent of the spin 
evolution. 

For /3 = we can make use of Eq. ^ and obtain a par- 
ticularly simple path integral representation for the spin 
evolution. Consider spin correlations with an initial den- 
sity matrix |cr)(cr| built from one of the two Kramers de- 
generate ground states of m = and cr — ±, having mo- 
menta Q = (7o.± = with G = i — + cP- . Using 
Eqs. (fTU]) and ^ for the forward and backward spin evo- 
lutions, we find an exact mapping of the spin correlations 
onto a superposition of equilibrium correlations of spin- 
less particles on a ring with a flux $ — iC (in units of 
quantum flux): 

Pa,*(t2i) = (e-'"^('=)e^''^('i))* . (11) 

Here again, 0(t) is a non-compact variable within (— cxd, oo) 
to be used within the path integral representation of the 
spinless problem. Note that the bath still couples to e^*^ 
hence we expect that (|11|) depends only on the noninte- 
ger part of $. For (5'a;(i)S'a:(0)) we obtain the following 
identity for an initial density matrix, 

Cl+{^ - isin2a(Fi,G(t)+F_i,G(i)) (12) 
+ cos^ |p_2g,g(0 + sin^ |P2-2G,G(i) • 

For C (t) the same result holds with all subscripts of 

Pa,Q reversing sign. For the correlations, on the other 
hand, we obtain 

C^+(<) = cos^ a + P-x-2G,G{t) sm^ a , (13) 

and for Cl_ (t) the same holds with Pi+2g,-g- Notice that 
the degeneracy point, a = -y/S, corresponding to fluxes 
$ = ±i represents a special case, and is not studied here. 

While the correlation function of the z-component of the 
spin, C^{t), obviously contains a constant non-decaying 
piece, the correlation function C^{t) contains only phase 
correlation functions Pa,<s> with a ^ 0. There is some ev- 
idence that these correlations decay in time. In particu- 
lar Pi ~ 1/t^ from the XY lattice model [16] and from 
small rj expansion [25] . Correlations with incommensurate 
a were studied in a related system of dissipative Josephson 
junctions [27], and found to decay algebraically. Further 
evidence is for large ry, as discussed below. To further 
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appreciate these correlations we have evaluated the path 
integrals in (jlip analytically for rj = 0, and surprisingly, 
we find P^2G ±g(^) ~ ^- a consequence, for 7y = the 
correlation function contains a piece which does not 
oscillate. As discussed below, though reduced, this part 
seems to survive for very weak dissipation, 77 <C 1, while 
it apparently decays algebraically for strong dissipation, 
ry > 1. 

Strong dissipation limit. — In the strong dissipa- 
tion limit, 77 ^ 1, we can describe the evolution of the 
phase through a Langevin equation, and an expansion in 
is possible [TTIHH]- In this limit, all correlation func- 
tions Pa,q with a ^ are found to decay algebraically. 
Large 77 perturbation theory yields that Pa,* ^ t^"" 
and the x-component of the spin also decays algebraically, 
while the z-component remains finite and does not de- 
cay. This result, however, holds only for up to times 
Int < 0(77) beyond which effects of renormalization of 77 
cannot be neglected. In a recent work [18 we have shown 
that in presence of a weak DC electric field there is a crit- 
ical 77c = l/27r such that 7/ > 77c flows to 7/c which would 
indicate Pa,* ^ i • In some sense the fluctuating spin 
corresponds to a time dependent flux, i.e. an electric field, 
though the correspondence is not precise. 

Weak dissipation. — The rather different behavior 
of Sx,y and Sz should already be manifest in the weak 
dissipation limit, where we can perform perturbation the- 
ory in the strength of the dissipation, 77. To do pertur- 
bation theory, we restrict ourselves to the case 5 = 1/2 
and /3 = 0, and use Abrikosov's pseudofermion method 
to represent each spinor ipma- of (l6|) by a pseudofermion 
operator, fma- In this language the ring Hamiltonian be- 
comes Hring = J2m,cr ^rn a Pm aJm a, while the interaction 
is expressed as 



(14) 



and standard field-theoretical methods can be used to 
evaluate physical quantities. A renormalization group 
analysis of the vertex function and the pseudofermions' 
self-energy reveals that, although ultraviolet logarithmic 
divergencies appear in both quantities, they cancel and 
the dissipation parameter 7/ is in leading order, neverthe- 
less, exactly marginal^ and the mass of the particle remains 
also unrenormalized [23_. 

In this perturbative regime, fingerprints of a non- 
exponential spin decay should appear in the susceptibility, 
X, which, in the absence of spin decay, should contain a 
Curie part. To compute x, we first express the impurity 
spin operator in terms of pseudofermions as 



EC« ft f 



(15) 



m,cr,m/,a"/ 



with the matrix elements simply determined from the wave 
functions dS]), as S'!^ ,,,^, ^, = mA^i\'^ ^n,a) ■ The lead- 
ing corrections to x are shown in Fig. [2] Although all 







Fig. 2; Zero and leading order (~ 77) corrections to the spin 
susceptibility. Continuous lines and wavy lines represent the 
pseudofermion and bosonic propagators, respectively, while the 
dots represent spin vertices. Logarithmic divergencies in the 
three diagrams above exactly cancel. 



corrections shown contain logarithmic ultraviolet singu- 
larities, remarkably, all these singularities exactly cancel, 
and one finally obtains just a finite renormalization of the 
perpendicular Curie susceptibility 



cos'^(a/2) 
4T 



The prefactor cos^(a/2) originates from the coefficient of 
P-2G,G{t) in Eq. ([T2|) . and accounts for g-factor renor- 
malization in the isolated ring. The correction ^ rj, on the 
other hand, represents the environment-induced renormal- 
ization of the X and y components of the spin (g- factor) . 
The above perturbative result and the survival of the Curie 
susceptibility indicates that the term P-2G.G{t) decays to 
a reduced but non-zero value for small rj. 

In contrast to the x and y components, the z component 
of the susceptibility, Xz , is found to remain unrenormalized 
by rj to leading order in the dissipation. These results 
imply that, for weak Ohmic dissipation, the only effect of 
dissipation is to slightly and anisotropically renormalize 
the 5-factor, but apart from that the spin behaves as a 
free spin, and does not decay. 

Case of /3 7^ 0. — So far we discussed only the case 
(3 = 0. We show now that the system with both a, /3 finite 
is equivalent to the Hamiltonian p^ . Assume a state \q) 
that is an eigenstate of Q\q) = q\q). This state generates 
a ladder of states, \m + q) = c""^\q), with integer m by 



Qe^"^'\q) = {m + q)e"^''\q) 



(16) 



Since f^^QT = —Q, a sequence of time reversed states is 
also generated by the time reversal operator, T: QT\m + 
q) = —{m + q)T\ni + q) . All these states are orthogonal 
since they correspond to different energy eigenvalues, and 
the environment couples the m and m ±1 states, exactly 
as for (3 = 0. The only difference is that Eo^{a,/3), which 
is not known analytically, changes the factor | — | vl-l-c? 
in E„if, Emi- Hence Eq. (fT4|) is a correct representation 
also of the /? 7^ caself] 



■^The matrix elements of the spin operators are nevertheless dif- 
ferent, changing e.g. the overall coefficient in Eq. II15I I. 
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Conclusions. — We derived the effective Hamiltonian 
of an electron confined to a ring within a 2-dimensional 
electron gas, in the presence of SO coupling, and subject 
to a dissipative environment. We have shown that the 
orbital motion of the particle decouples from the spin evo- 
lution, and correspondingly, spin decay has a geometric 
character [T]. For an Ohmic environment, we mapped the 
spin relaxation problem to that of a spinless particle on a 
ring pierced by a magnetic flux [Eqs. (|12I13|) ] . We find that 
the z component of the ground state spin is not affected 
by dissipation. The x and y in-plane spin components are, 
on the other hand, reduced by dissipation, but we find no 
dephasing for spin | and weak dissipation. However, these 
components seem to dephase at large dissipation. 

We should remark that these latter results are based on 
the assumption of Ohmic dissipation. The situation may, 
however, change for subohmic dissipation or Xjuf* noise, 
present in many systems. In this case, the decoupling of 
the spin and orbital motion and thus Eqs. (I12I13P remain 
valid, however, for subohmic dissipation 77 is a relevant 
perturbation, and even a small dissipation could possibly 
lead to the decay of the x and y spin components. This 
possibility, however, needs to be further explored. 
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